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Abstract 

Generalization of Gaussian trial wave functions in quantum molecular dy- 
namics models is introduced, which allows for long-range correlations char- 
acteristic for composite nuclear fragments. We demonstrate a significant im- 
provement in the description of light fragments with correlations. Utilizing 
either type of Gaussian wave functions, with or without correlations, however, 
we find that we cannot describe fragment formation in a dynamic situation. 
Composite fragments are only produced in simulations if they are present as 
clusters in the substructure of original nuclei. The difficulty is traced to the 
derealization of wave functions during emission. Composite fragments are 

produced abundantly in the Gaussian molecular dynamics in the limit Ti — > 0. 
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I. INTRODUCTION 



The theoretical description of heavy-ion collisions at intermediate beam energies, 10 MeV 
~ E/A ~ 150 MeV, is still in an unsatisfactory state. Different factors contribute to that 
situation. Thus, for one, the nucleon excitation energies are low and in that comparable to 
nucleon localization energies, indicating a likely importance of the quantal effects. With the 
change in the beam energy in the discussed range, the dynamics changes from that dom- 
inated by the mean field to that dominated by collisions (as evidenced in the appearance 
of the flow balance-energy). As excitation energies grow, they begin to exceed the average 
binding energies per nucleon and, within the mentioned range of E/A, a massive production 
of intermediate- mass fragments (IMF) takes place 0. The production, in terms of IMF 
multiplicity or total mass that IMF carry, maximizes at E/A ~ 75 MeV. The description of 
the intermediate and light fragment production is beyond the capability of common single- 
particle models of collisions 0. The single-particle models with fluctuating forces |||[] can 
describe fragment production, but miss the shell effects and the discreteness of the mass and 
charge numbers. The involved limitation is recognized once one realizes that, in the very 
central Au + Au collisions at 100 MeV/nucleon, the probability for a proton to come out 
from the reaction as a constituent of an a particle is close || to 50%. Within the Boltzmann- 
Langevin model with the fluctuating forces [0,0, the a particle plays no distinguished role. 
Statistical models PJ^] account for the shell effects, but miss the reaction dynamics. The im- 
portance of the dynamics is seen, in particular, in the fact that the collective outward flow 
energy in the reactions is comparable to the thermal energy. 

Given the above situation, the quantum molecular model proposed in recent years for 
the reactions [§ was met with quite some expectations. In this model, the quantal wave 
function for a reacting system is represented as a product of Gaussian wave-packets for 
individual nucleons. The packets have dynamic centers, phases, and widths. The parameters 
obey equations of motion following from a variational principle. The model || accounts 
for shell effects and has been shown to describe the evaporation of individual nucleons 
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from excited nuclei 0. As nucleons leave a nucleus, their wave-packets become completely 
delocalized. The packets are used in this model in favor of the Hartree-Fock wave-functions, 
because with the packets one expects to describe fluctuations. Note that initialization of 
the classical Vlasov equation with a set of 5-functions leads to the same results as the exact 
equations of motion. The set of wave-packets is the closest approximation to the set of 
5-functions that one can get quantally. Other molecular models in use might be considered 
amended or simplified versions of |§. Both in the model |§ (termed FMD) and in |To|| 
(termed AMD), the overall wave function is antisymmetrized. However, in |10| the dynamics 



of the wave-packet width is suppressed. In |TT| (EQMD) the width is dynamic, but the effects 
of antisymmetrization are accounted for only approximately, using a Pauli pseudopotential. 
In [12| and |13| (QMD) neither the widths are dynamic nor the antisymmetrization is carried 



out explicitly. These models are the most classical within the class. In the models AMD, 
EQMD, and QMD, collisions between wave packets take place, on top of the wave-packet 
dynamics obtained or attributable to the variational principle for wave functions. 

Within the molecular models, fragment production has been studied, in quite some 
detail, in the most classical of the class, the QMD model |T3,r5|. Within the FMD model, 



the fragment emission has been observed in the calculations of reactions involving light 
nuclei |§ . In the FMD model a specific problem arises concerning the emission of fragments. 
Inside a fragment, the constituent nucleons are localized in the relative separations. When 
the wave function is a product of single-particle wave functions, this implies a localization 
of the fragment center of mass. While low-energy nucleons escaping from a nucleus become 
delocalized, getting rid of the kinetic energy associated with their localization ||, this cannot 
be the case for the fragments. The cm. localization energies e.g. in the range of A = (2- 
4) nuclei, are within AE ~ (10 — 30) MeV and, given that the temperatures in excited 
nuclei in reactions could be as low as (5-10) MeV, they could result in significant thermal 
penalty-factors for emission, e _AE//T . 

Generally, within FMD and EQMD the kinetic energy of the localization of the center of 
mass is not conserved and may be transformed into the internal energy, in a translationally 



invariant situation. While this need not be a problem for large fragments with a small 
localization energy compared to the fragment total energy, for light fragments it can mean 
that their excitation energies cannot be determined. 

In this paper, we investigate fragment production within a model of the FMD type. 
Given the problems associated with the localization, mentioned above, we consider the wave 
functions, of a Gaussian form, that are products of wave functions for individual nucleons 
and also such that allow for the correlations between nucleons within fragments, with a 
decoupling of the fragment center of mass motion. At present, effects of antisymmetrization 
are accounted for approximately only using a Pauli potential. 

The outline of the paper is as follows. Section |I| describes the trial wave functions for 
the Schwinger variational principle. The equations of motion following from that principle 
are discussed in Sec. |T|. Section [IV| discusses a choice of the hamiltonian which permits 
an analytic calculation of the expectation values in the equations. Our results on fragment 
dynamics and production are reported in Sec. [V] and the conclusions are given in Sec. [VT|. 

II. TRIAL WAVE FUNCTION 

The trial wave-function for the Schwinger variational principle is taken in the form 

(aji, ...,x N \V) = Afexp (</>Oi, • • • , en)) \ , (1) 

where x is a normalized spin-isospin wave function. The argument in the unnormalized 
exponential wave function is 

<f>(xi, . . . , x N ) = -Aij fa - ri) ■ (xj - Vj) + ipi Xi. (2) 

Here, £Cj denotes the position vector of particle i. The repeated particle indices indicate sum- 
mation. In absence of antisymmetrization, from ([[]) and (|2|) one finds that the normalization 
constant is equal to Af = (det(2Rev4)/7r Ar ) 3 / 4 , where N is the particle number. 

The parameters of the spatial wave function, which depend on time, include the elements 
of the complex symmetric 3iV x 3iV matrix A with a positive definite real part, and further 
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the components of the 2N real vectors {r^} and {Pj}. Below, we shall use q u to indicate any 
member of a set of the dynamic parameters. A convenient mapping of the set follows using 
Qi(i-i)/2+j — ReAjj and qN(N+i)/2+i(i-i) /2+j = Im A^, where in the both equations j < i, and 
further using q N{N+1)+3 ^i )+a = (n) a and gAr(Ar + 4) + 3(;_i) +a = (pA a . The time-dependent 
parameters relate to different expectation values with 

(xi) = ri, (-iVi)=Pi, (3) 

((«.-)« («i) 6 > = (r0« (r-i) 6 + ^ (ReA)r; (4) 

and 

(i{Vi)ai(Vj)b) = (Pi)« (Pj)fc + *«6 (Re A)- 1 A)jj , (5) 

where a and b are indices for carthesian coordinates. Given the parameters r*j and 

Pj are referred to, further, as the centroid and momentum of a particle i, respectively, and 
A is referred to as the width matrix. 

The form (|I|) includes the special case of the width diagonal in particle indices, Aij = 
Sij Ai. The wave function in that case reduces to a product of single-nucleon wave functions 
such as utilized in the FMD or EQMD calculations. To illustrate the advantage that ([!]) 
with (0) offers, let us consider a deuteron. The parametrization allows for the wave function 
of the form 

(xi, x N \^) = A/"exp ( — A cm {xt + x 2 ) 2 /4 - A rel { Xl - x 2 ) 2 ) x , (6) 

where A cm = 2 (An + A12), A re i = (An — Ai 2 )/2, and An = A 2 2- The magnitude of 
A re i can be adjusted to reproduce the r.m.s. radius of the deuteron and A cm in the above 
may take on arbitrarily low values corresponding to a delocalized deuteron, as expected in 
the emission in reactions. FMD and EQMD parametrizations more standard than ours, 
with vanishing off-diagonal terms in A, on the other hand, permit only strongly localized 
deuterons, as A cm = 4A re i for A\ 2 = 0. 
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III. EQUATIONS OF MOTION 



The wave function (|IJ) depends on time only indirectly through the parameters q^. In 
this section, we shall obtain equations of motion (eom) that govern the behavior of these 
parameters. 

On writing the time-dependent variational principle in the form 5 J tl 2 C{q^^ q^)dt = 0, 
the eom for g M follow (formally) as Lagrange- Euler equations, — J^- = 0. Using 

C = (^\ih^r — H\ty) for the Lagrange function, one obtains, see [(J, 



•Avfi q^i 







(H) 



(7) 



where the matrix A, multiplying the time derivatives, is skew symmetric. This matrix is 
related to the overlap of the derivatives of the wave function with respect to the parameters: 

d 
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(8) 



Given the wave function of the form (|1|), one finds for A 
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dq. 



(9) 



where the square brackets stand for [O] = J HfLi dxi exp (4>*)0 exp (0). From (|), one can 
see that A does not couple parameters describing spin-isospin degrees, with those in the 
spatial wave function. For an interaction diagonal in spin and isospin, considered below, 
this implies that the spin-isospin wave-function does not depend on time. Correspondingly, 
the spin-isospin wave functions will be largely disregarded in the further discussion. 
The explicit expression for A, obtained using (|2]), is 

" = 4 (Re A) - (Re ^— + Wu W, " {v " ^ ' (10) 

Since A is independent of r*j and p { , the eom for particle centroids and momenta take on the 
form of the Hamilton's equations with the expectation value of the Hamilton operator in 
these equations playing the role of a classical Hamiltonian. As far as the width is concerned, 
for a practical solution of the eom, it is necessary to invert the ^4-matrix. When represented 
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in the space of the real and imaginary elements of A, the matrix A acquires a 2 x 2 block 
structure with the blocks on the diagonal vanishing and those off-diagonal (coupling the real 
and imaginary parts of A) given by ±B. Fortunately, the quadratic matrix B, of a dimension 
iV(iV + l)/2, can be easily inverted analytically with the result 

^n(n-i) . i(i-i) ■ = 77 (Re^niRe^mj + Re A n ~ Re A m i) , (11) 

v 2 ' +m, 2 +J 3 

where 1 < m < n < N, 1 < j < i < N. The equations of motion for the parameters - 
centroids, momenta, and matrix elements of A - take then on the form 

9(H) . d(H) 

dpi dr { 

j t Re Aij = E n B 5k +j « +m d i^ Anm > 1 Im A = ~ B ~ l • (13) 



The indices in the second equation in (|13|) should be handled in the same manner as in the 
first of the equations where they are written out explicitly In evaluating the derivatives 
in fljjD, the elements ReA nm and ReA mn , and lmA nm and lmA mn , respectively, should be 
treated as identical. 

Given (|H|), the contributions from the kinetic energy to the eom ( |T3| ) for the width 
matrix may be simplified into 

B-.m^w, B-^m- ™ ReA *. (14) 

a 1m A m a Re A m 

For the free particles, eom may be then represented as A = —2ihA 2 /m, with the solution 
A(t) = Aq (1 + 2ih(t — t )A /m)~ 1 . This means that, for free particles, a width matrix 
diagonal initially in certain directions in the iV-particle space will stay diagonal in these 
directions. Contributions to the kinetic energy associated with different directions, propor- 
tional to |v4(t)| 2 /Re A(t) = |j4 | 2 /Re A , stay constant as a function of time. For interacting 
particles, clearly, interaction terms appear on the r.h.s. of the eom for the width. These 
terms can, generally, cause violations of the conservation law for the cm. kinetic energy, 
when the width matrix is constrained to the diagonal form. 



IV. HAMILTON! AN 



The eom (12) and (|i~3|) express the time derivatives of the parameters q a in the wave 
function as the linear combinations of the derivatives of the expectation value of the Hamil- 
tonian, with respect to q a - We shall now calculate the Hamiltonian expectation value in 
terms of the parameters. 

The expectation value of the kinetic energy, (T) = -^£i(V7), is directly obtained 
from (H). To make the calculation of the expectation value of the potential analytic, we 
choose the internucleon potential of such a form as in ||, 

V = X>*E(«'* + ( 1 - v) k )pf) exp ( - A, 2 ( Xl - x s f) , (15) 

k i<j 

where P^ is the Majorana operator exchanging spatial coordinates of the particles % and j. 
The sum over k extends over the repulsive and attractive contributions to the potential. 
The values for Uk, \k and Wk are listed in ||. With (|15|) , we exclude the Coulomb potential 
as in |J, but then we only consider the light to medium nuclei. The expectation value of 
the potential fll5|) with respect to the wave function ([j]) is then 

(V) = E^E^ 2 [exp {-K 2 { Xl - X] f) (w k + (1 - w k )P%) } (16) 

k i<j 



where the meaning of the square brackets is such as below Eq. (^). On defining the matri- 
ces A* J and Ap by the relations, respectively, A^ m r„ ■ r m = (r* — rj) 2 and A p m = A P ! n ) p( m ) 
where P exchanges indices i and j, we obtain 



exp(-A fc fa - Xj) )0 = -372 exp(y# ) . (17) 



(det(A + A°* + Afc 2 A 

The argument of the exponential on the r.h.s. of ( |T7|) for the two cases of the operator O is 



and 



jj,k (A + A p * + A fc 2 A lj ) n m ^ P ^ P (A + A p *) nm ^ P ^p 

<PO=P — a In Qm ' A S n S m 1 l iy J 



with s£ = r n - r P(n) and g£ = p n - p P{n) - i(A - A p *) nm s?. 

At first glance, the calculation of (V) for non-diagonal width-matrix A scales with the 
number iV of particles as N 5 ; for each of N(N — 1)/2 particle pairs, one has to invert a. N xN 
symmetric matrix. Fortunately, for each of the particle pairs the matrix that needs to be 
inverted differs from the matrix 2 Re A only in 4N — 6 elements. Thus, for each of these 
pairs one can use information from inverting the matrix Re A, and the calculation of (V) 
scales then only with the particle number as N 4 . 

Generally, one would want the trial wave function to be antisymmetrized. However, 
the width matrix A nondiagonal in particle indices introduces such a number of parameters, 
that the antisymmetrization ceases to be feasible for the particle number larger than ~ 16. 
Thus, we resort to the Pauli potential to simulate the effects of antisymmetrization. The 
Pauli potential acts only between particles with the same spin and isospin, and it is chosen 
proportional to the Majorana operator, Vp au u = u$ P^ . The latter is motivated by the 
fact that the expectation value of P M for the product wave function is proportional to the 
square of the scalar product of single-particle wave functions, {P-f) = where 
ipij are the single-particle wave functions of the two particles. When one wave function 
approaches another, the system reacts with a repulsive force. The Pauli potential is added 
as the third component to fll5f) (1U3 = 0, A^ 1 = 0), and the value of U3 is adjusted to best 



fit the properties of ground-state nuclei. Unfortunately, such w 3 depends quite significantly 
on the mass. Using parameter set Bl from || (originally from | |16|| ) in (|15D, we obtain 



e.g. optimal w 3 = 70 MeV for A = 12 and % = 200 MeV for A = 80. (This would have been 
likely partly alleviated if we included the Coulomb potential.) 

As an example, in the obtained ground state of 12 C three four-fold degenerate centroids 
position themselves at the corners of an equilateral triangle in configuration space. The av- 
erage associated particle momenta are zero. For the diagonal width matrix and the Bl 
interaction, the minimum internal energy (defined as the difference between the total en- 
ergy and the energy of the center of mass motion) of E int = —87.1 MeV is obtained for 
Re An = 0.33 fm~ 2 which gives R rm s = 2.32 fm. The obtained 12 C nucleus is stable against 
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break up into three a particles which have a ground state energy of Ei nt = —28.1 MeV. (Note 
that the 12 C energy does not contain the Coulomb energy estimated at Eq = 11.4 MeV using 
the formula (61) in |J; the value of R rm s includes the spatial extent of a proton ||.) For the 
width matrix with off-diagonal elements, a lower minimum of Ei nt = —95.4 MeV is obtained 
for Re An = 0.30 fm~ 2 , and Re A^- = —0.038 fm~ 2 , if j is within the same cluster as i or has 
the same spin and isospin directions as i, and Re A%j = 0, otherwise. There is a freedom, 
in the latter case, of adding a constant to all elements of A that only changes constraints on 
the cm. motion. 

For even-even nuclei of mass larger than carbon, our Pauli potential favors differences in 
average single-particle momenta over differences in centroids in the ground state. The cen- 
troids for these nuclei become identical while the momenta get distributed in momentum 
space in the groups of four. For example, in the ground state of 16 the momenta are placed 
at the corners of a tetrahedron in momentum space. 

V. SOLUTION OF THE EQUATIONS OF MOTION 

We first discuss differences in the dynamics of isolated light fragments, for correlated and 
uncorrelated trial functions. We then investigate fragment production within the present 
dynamic description, when a compressed and excited nuclear system expands and when 
nuclei collide at low energies. Some of our results are quite unexpected. 

Since the differences in the dynamics for correlated and uncorrelated trial functions are 
expected to be the largest for the lightest of fragments, we investigate the dynamics of 
an isolated deuteron and of an isolated a particle, illustrated in Figs. 1 and 2, respectively. 
As the deuteron is unbound for all the interactions listed in 0, we use the Volkov 1 (VI) |17j 
with u a = —104.5 MeV in the deuteron case. Fig. la displays the evolution of total E to t, 
kinetic E^in, potential E pot , and total internal Ei nt energies, for a deuteron initialized in 
the state of the lowest internal energy. Either the dynamics for a correlated (labeled c) or 
uncorrelated (labeled u) trial-function is followed. In the case of the correlated dynamics, 
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both the total and the internal energies remain constant as a function of time. However, 
in the case of the uncorrelated dynamics, the internal energy increases, at a cost of the center- 
of-mass energy. At t ~ 12 fm/c the deuteron becomes, in effect, unbound and it remains 
so thereafter. The different evolution is associated with the behavior of matrix elements, 
as illustrated in Fig. lb. In the absence of correlations between particles, i.e. A 12 = in 
Eq. (fl), the derealization of the deuteron center of mass couples to the derealization of the 
internal state. In particular, a reduction in Re A cm requires a reduction in Re A\ and this 
implies a reduction in Re A re \. By contrast, in the correlated dynamics, the width for the 
center of mass, related to Re A cm = Re Aij, behaves like the width of the Gaussian wave- 
packet for a particle with twice the nucleon mass. While the element A cm decreases with 
time, it does so because the off-diagonal element A 12 of the matrix increases in magnitude. 
The element A reh cf. below (||), stays constant. 

As we have already demonstrated, for the uncorrelated dynamics an unphysical exchange 
of energy occurs between the internal degrees of freedom and the center of mass motion. 
When the magnitude of the internal energy exceeds the energy of localization of the center 
of mass, the coupling may cause unphysical oscillations for the ground state. This is shown 
for the a particle in Fig. 2. In the case of a correlated wave function, the width for the 
center of mass behaves as the width of a Gaussian packet for a particle with four times the 
nucleon mass. The internal part of the wave function does not change with time. In the case 
of an uncorrelated wave function initialized in the lowest state of internal energy, the matrix 
element and, correspondingly, the width for the center of mass oscillate. An exchange of 
energy, back and forth, occurs between the internal and center-of-mass degrees of freedom. 
If the uncorrelated wave function is initialized in the lowest state of total energy, the internal 
and center-of-mass energies stay constant. The internal wave function does not vary with 
time, but nor does vary the center-of-mass wave-function, with the center of mass never 
getting delocalized. Needless to say that in that case the internal energy is higher than in 
the ground state. The examples in Figs. 1 and 2 show the benefits of using the correlated 
over the uncorrelated wave functions in the fragment description. 
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Turning now to fragment production, we start out by exploring the situation where a 
highly excited system formed in the central heavy-ion collisions decays into vacuum. For such 
a system we expect a reduced importance of the antisymmetrization that is missing from 
our equations. We assume that internal degrees of freedom of the system are, generally, 
equilibrated, allowing only for a variable strength of the radial flow. We investigate the 
dynamics with a variable width matrix, either restricted or not to the diagonal form and, 
further, the dynamics with a static width, i.e. classical. More significant and surprising 
differences are found between the quantal and classical dynamics than between the dynamics 
with different variable width matrices. 

To simulate the excited system (A = 20 — 80), we distribute centroids randomly within 
a spatial volume of radius R and in the momentum according to a finite-temperature Fermi 
distribution (T = 5 — 12 MeV). To account for the flow, we add to the momentum of each 
particle a component proportional to the position vector relative to the overall center of 
mass. The proportionality constant determines the amount of flow energy in the initial 
state (E col = — 25 MeV/nucleon). The width matrix is initialized as a real multiple of the 
unit matrix. Spatial representation of one of such initial states is given in Fig. 3a. In the 
particular case A = 80 and R ~ 5 fm. The rms radii of individual packets are equal to 1.9 fm. 

In the case of a dynamic matrix, whether or not restricted to a diagonal form, the ex- 
cited system, initialized as above, emits a number of single nucleons in the course of time. 
The number of emitted nucleons generally increases with the energy of the system. How- 
ever, at no particular energy, for the studied A = 20 — 80 systems, any emission of IMF or 
even of a particles is observed. This appears to be true irrespective of how we divide the 
excitation energy into collective and thermal. Examples of the late-stage distributions of 
centroids in space are given in Figs. 3(c) and 3(d) for the initial net excitation energy of 14 
and 26 MeV/nucleon, respectively. For the particular initial states, the centroid distribu- 
tions at the respective time are not distinguishable by eye between the evolutions with and 
without correlations. We found that frequently to be the case for the systems initialized in 
the manner discussed above. While some centroids appear close to each other in Figs. 3(c) 
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and 3(d), at a distance from the main residue, the respective wave-packet widths are so 
large that the packets will separate from each other, eventually. At t — 300 fm/c, for the 
displayed systems, all relative energies of the emitted particles are positive. 

To check whether the lack of fragment production in the correlated dynamics might 
be associated with the initialization of the dynamics in an uncorrelated state, we have 
carried out tests using different initializations. Thus, we have added to the hamiltonian 
an oscillator term V = v osc to keep an excited system from expanding into vacuum, 

while allowing the matrix elements of A to thermalize, eventually removing this additional 
potential. Figure ^ shows the late stage of an A = 80 system initialized using a narrow 
oscillator potential. During the time of 100 fm/c within the oscillator potential, the rms 
values of off-diagonal elements of the matrix A had saturated, ((Re A^j) 2 ) 1 / 2 ~ 0.7 /(A — 
1); for the chosen v osc the system hardly expanded in the potential, compared to t = 0. 
On removing the constraining oscillator potential, we added collective components to particle 
momenta, proportional to the distance from the center. The system, subsequently, emitted 
a number of single nucleons and a deexcited residue formed at the center, consisting of 37 
nucleons at the time shown in Fig. |j. In tests we changed the extension of the constraining 
oscillator potential for the excited nucleus, the initial temperature, and the magnitude of 
collective energy. Further, we initialized excited nuclei without any initial constraining 
potential, just assigning random gaussian values to the off-diagonal and diagonal elements 
of A, using the thermalized values from the oscillator potential as a guidance. Consistently, 
in all tests, the released nuclei emitted, in the course of their evolution, a number of single 
nuclei but never any IMF or even an a particle. Thus, the particular feature does not depend 
on the off-diagonal terms of A being zero or finite in the initial state of an excited nucleus. 

Clearly, in the past, the production of IMFs and light clusters has been observed in the 
QMD calculations. The QMD limit corresponds to taking h —>■ 0, or to suppressing the 
width dynamics in our equations. Indeed, when taking h reduced by a factor of (5-10), 
or the width dynamics slowed down by a such a factor, we begin to observe the cluster 
production. Figure 3(b) shows centroids for a system evolved from the initial state shown 
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in Fig. 3(a), using a static width, for the same available energy (given the frozen width) 
as the system in Fig. 3(c). At a time t = 300 fm/c in Fig. 3(b), two IMFs as well as 
two dinucleons are seen and they will remain stable. Additional small clusters have been 
emitted before t = 300 fm/c and left the displayed spatial region. As the reason for the lack 
of cluster production in the calculations with a dynamic width matrix (whether diagonal or 
not) emerges, in our tests, the spreading of wave functions towards the emission time. After 
packets get delocalized, the interaction is not capable to contract them back into fragments. 

Our results on the cluster production for the dynamic width matrix may seem in con- 
tradiction to the FMD results with even multifragmentation events reported in nuclear 
collisions. However, a scrupulous examination indicates that the clusters, seen in the final 
states of the FMD calculations in ||, were not formed during reactions, but were present 
in the initial states and survived reactions. When investigating that particular issue, we 
simulated symmetric collisions of nuclei with different initial structures. Thus e.g. for the 
potential ( |i~5| ) and our Pauli potential, the centroids within the ground state of l2 C form 
three a-type clusters of four nucleons each. On the other hand, within 40 Ca the nucleon 
centroids situate themselves at the overall center of mass position; the widths for different 
nucleons take on different values. In the true ground state for our potential, the momenta 
form groups of four in 40 Ca, with the four particles being two protons and two neutrons 
with different spin directions. To illustrate, though, our point on reactions, we shall displace 
slightly the momenta from the identical values in 40 C a in the initial state of a reaction, 
making sure that change in the overall binding energy is negligible. Notably, inclusion of 
any kind of spin-isospin dependent interaction would break the momentum sub-clusters in 
nuclei of A < 40, anyway. Then, when considering the 12 C + 12 C reaction, we shall deal 
with six sub-clusters in initial state, and with none in the m Ca + 40 Ca case. 

Figure [5] shows the initial and late states of an exemplary 12 C + 12 C reaction at the beam 
energy of 29 MeV/nucleon. While one nucleus got highly excited in the reaction, the other 
has fragmented into three a particles. Each of the a particles is excited at the displayed 
time, however remaining below threshold for particle emission. 
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Figure ^ shows three stages of a 35 MeV/nucleon m Ca + 40 Ca reaction with no sub- 
clusters in the initial state. At the first of the times shown, the wave functions of the two 
nuclei just started to overlap. At the second of the times, a transient residue that formed is 
maximally spread out. At the third of the times, the outcome of the collision is, essentially, 
determined and only a residue and some single nucleons are seen. No clusters are emitted 
from this collision. 

The examples presented in Figs. ||] and ^ are typical ones: for none of the studied initial 
states, we have observed, for the dynamic width, the emission of IMF's or of a particles that 
were not present in the substructure of an initial state. This has been the case whether we 
included or excluded the correlations. Again, as a reason for the absence of cluster formation 
for the dynamic width, we find that, at times when the system reaches a density when clusters 
are expected to form, the width of the wave packets has grown so large that the interaction 
between different wave packets is too weak to force the width of these packets to shrink, 
and their centroids to approach each other enough to form a nucleus (see also [[HJ). When 
we suppress the dynamics of the matrix A, new clusters form and emerge from the reaction 
region. 

The values of the widths are not necessarily unphysically large, as nucleons are expected 
to get delocalized with the reaction progress. In reality, though, when broad wave packets of 
different nucleons overlap in space, the interaction between the nucleons should be generally 
able to clump two or more nucleons, within a distance comparable to the interaction range, 
into a cluster. The wave function in terms of broadened Gaussians (whether or not diagonal 
in the nucleon coordinates) does not allow for that. 

VI. SUMMARY AND CONCLUSIONS 

We have investigated the time evolution of nuclear systems, in terms of correlated and un- 
corrected Gaussian wave functions, following from the time variational principle. As an in- 
teraction in the Hamiltonian, we utilized a Volkov-type potential. Upon determining that 
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an antisymmetrization of the correlated wave function would not be feasible beyond rela- 
tively small systems, we employed a Pauli type potential to simulate the effects of Pauli 
principle. 

For uncorrelated trial wave functions, the internal state cannot be localized without 
localizing the center of mass. Thermal estimates indicate that this could suppress fragment 
emission. In dynamics, an unphysical energy exchange occurs between the center of mass 
and the intrinsic motions. These deficiencies are absent when using the correlated wave 
functions. We have explicitly demonstrated an improvement, in terms of correlated wave 
functions, in the description of deuterons and alphas. 

Contrary to expectations, the inclusion of correlations has not improved the description 
of cluster production in the explicit simulations of heavy-ion reactions or in situations char- 
acteristic for the reactions. Either in correlated or uncorrelated Gaussian dynamics, clusters 
are only produced when they are present in the substructure of an initial state. This appears 
to be true, in the correlated dynamics, whether or not the initial state is correlated. The ab- 
sence of new clusters is associated with a large spreading of the wave function at the time 
when a reacting system expands and the new clusters are expected to form. The spreading 
in relative coordinates is present even in the correlated wave function, despite of the fact 
that, for a correlated wave function, the relative spreading may evolve independently from 
the spreading for the center of mass for any of the potential fragments. The interaction is 
too weak to pull back the wave function to a size appropriate for a fragment. In classical 
dynamics fragments are produced. 

While there is nothing unphysical in the spreading of the wave function as such in the 
simulations, in reality the interaction would be capable of creating correlations in the wave 
function over distances of the interaction range. E.g., when considering the two-nucleon 
correlation function, with time the function should develop into a spreading long-range 
component, weighted with a certain probability, and a more stable short-range component 
corresponding to the forming fragments. The specific parametrization of the wave function, 
once the wave function has spread, does not allow for the development of the short-range 
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component. 

On a general level our results show that the dynamics is important in fragment produc- 
tion. The spread of the wave packets at emission indicates that an emitting system can 
gain the packet derealization energy. On the side, within AMD the width dynamics is sup- 
pressed and an ad hoc term is added to the Hamiltonian to account for the derealization 
at emission. Within QMD, the location energy is disregarded at every stage, including the 
initial ground-state nuclei and the final state. Practically, the procedure in AMD amounts 
to a change in the effective interaction, compared to QMD. Within either the approaches, 
the fixed width of a wave packet extends the range of the two-particle interaction and may, 
in fact, act to suppress fragment formation. 

Recent results from FMD with short-range correlations [19j show some fragments in the fi- 



nal state of reactions, in contrast to the present calculations (or ||18|| ). This can be attributed 
to a lesser spreading of the wave packets towards emission time, due to the combination of 
the effects of short-range correlations and the assumed spherical shape of packets. Given 
that this lesser spreading is unrelated to the physics of fragment production, we believe that 
it is not a solution to the the dilemma at hand. In particular, we expect too low fragment 
yields, from the calculations with short-range correlations, compared to experiment. 

We conclude that, in a successful description of fragment production, the wave function 
or, more generally, the model density matrix must a have a flexibility to change over dis- 
tances comparable to the interaction range, at a time when the fragments are formed. One 
possibility within the dynamics is to keep the nucleon wave-packet width static and com- 
parable to the interaction range, and to account for a quantal spreading with time through 
a stochastic decision process This can be nominally derived through the reduction 

of the wave function space to the space spanned by the wave packets of constant width, 
which results in a residual force associated with the kinetic-energy part of the hamiltonian. 
The present authors have been exploring a replacing of the single wave packet for every 
nucleon by a superposition of packets. Unfortunately, this looses the inherent simplicity 
of FMD. The high hopes associated with the FMD approach cannot be directly realized as 
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far as fragmentation is concerned. 
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FIGURES 

FIG. 1. Evolution for the internal ground state of a free deuteron. (a) The dashed, solid, 
dash-dotted, and dotted lines show, respectively, the evolution of the kinetic, total, internal, and 
potential energies in correlated (c) and and uncorrelated (u) dynamics, (b) Evolution of the 
elements of the width matrix. The dashed lines shows the evolution of the element in the case 
of an uncorrelated wave function. The dash-dotted, solid, and dotted lines show, respectively, 
the evolution of the internal, center-of-mass, and diagonal elements in the case of a correlated wave 
function, cf. Eq. (||) and text below. 

FIG. 2. Evolution of the elements of the width matrix A for an a particle initialized in its 
ground state. The solid, dash-dotted, and dotted lines show, respectively, the evolution of the 
center-of-mass, relative between-two-nucleons, and diagonal elements, in the case of a correlated 
wave function, cf. Eqs. (^) and (Q). The long-dashed line shows the evolution of the element 
in the case of an uncorrelated wave function initialized in the lowest state of internal energy. 
The short-dashed line shows the evolution of the element in the case of an uncorrelated function 
initialized in the lowest state of total energy. 

FIG. 3. Particle centroids (dots) in the configuration space for A = 80, in the initial state (a) 
and at t = 200 fm/c for the static wave-packet width (b), and for the dynamic width matrix (c) 
and (d). In (c), the initial flow energy is lowered, compared to (b) and (d), by the initial energy 
content in the localization of the wave packets. The circles indicate rms radii of the most and least 
localized Gaussians. The radii of the packets in the case (b), and the radii of the most localized 
packets in the case (c), are comparable to the radii of the dots representing centroids. The axes 
show distances in fm. 
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FIG. 4. Late stage (t = 300 fm/c) of an A = 80 system initialized at a temperature T = 8 MeV 
and E co i = in a narrow oscillator potential V = v osc J2i x h where v osc = 1 MeV/fm 2 . At a 
time t = 100 fm/c, the oscillator potential was removed, and particles given collective outward 
velocities corresponding to net E co \ = 7 MeV/nucleon. The centroids for emitted and bound 
particles (negative and positive removal energies, respectively) are indicated with small open and 
small filled circles, respectively. The large circles indicate widths, from the diagonal elements of A, 
for most localized and delocalized bound and emitted particles. 

FIG. 5. Initial (left side) and late (t = 300 fm; right side) states of a 12 C + 12 C reaction 
at 29 MeV/nucleon. The beam axis is directed vertically, dots represent the centroids of nucleon 
packets, and circles show the packet rms radii. Edges of the spatial boxes show distances along the 
carthesian axes in fm. 

FIG. 6. Three stages of a 35 MeV/nucleon 40 Ca + 40 Ca reaction in the configuration space. 
From left to right, the boxes show the reaction, respectively, at t = 20 fm/c, t = 100 fm/c, and 
t = 300 fm/c. The dots indicate centroids, while the circles indicate the rms radii of the least and 
most localized Gaussians. Box edges show the distances along the carthesian axes in fm. 
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